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\(-} , We study the non-forward matrix elements of twist-2 QCD light-ray operators and 

their representations in terms of skewed and double distributions, considering the 
S— : ■ pion as well as the nucleon. We point out the importance of explicitly including all 

twist-2 structures in the double distribution representation, which naturally leads 
to a "two-component" structure of the skewed distribution, with different contri- 
r S^ butions in the regions \X\ > £/2 and \X\ < £/2. We compute the skewed and 

double quark distributions in the pion at a low normalization point in the effective 
C^h! chiral theory based on the instanton vacuum. Also, we derive the crossing relations 

expressing the skewed quark distribution in the pion through the distribution am- 
plitude for two-pion production. Measurement of the latter in two-pion production 
in 7*7 and 7*iV reactions could provide direct information about the skewed as well 
rS ' as the usual quark/antiquark-distribution in the pion. 
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1 Introduction 

In so-called light-cone dominated hard scattering processes the non-perturbative infor- 
mation entering the scattering amplitude is contained in matrix elements of certain QCD 
light-ray operators between hadronic states. A well-known example is inclusive deep- 
inelastic scattering, where in the asymptotic regime the cross section is determined by 
forward (diagonal) matrix elements of twist-2 light-ray operators in the target state, 
which have an interpretation as parton distributions. More recently, factorization has 
been proven for a large class of exclusive processes, namely deeply-virtual Compton scat- 
tering (DVCS) and hard meson production [|], ||, ||, |], ||, §, [7], ||. The amplitudes for these 
processes involve non-forward (more generally, non-diagonal) matrix elements of light- 
ray operators between incoming and outgoing hadron states, which can be represented 
by generalized parton distributions. Such matrix elements had earlier been introduced in 
the description of Z° photoproduction at small x and in the context of the non-local 
light-cone expansion [l(| . 



Due to the presence of a non-zero momentum transfer, non-forward matrix elements 
of light-ray operators possess a much richer structure than the forward ones defining 
the familiar parton distributions. In the matrix element of a generic light-ray operator, 
(p — r/2\if(—z/2) . . . if(z/2)\p + r/2), with z 2 = 0, both the momentum transfer, r, and 
the average of initial and final momenta, p, in general have non-zero longitudinal ( "plus" ) 
component with respect to the light-cone direction defined by z. In a partonic language, 
one may express the momentum of the "active" parton in terms of any linear combination 
oip-z and r-z. Two approaches have been proposed. One can analyze the matrix element 
assuming proportionality r-z = £,p-z, where the value of £ is determined by the kinematics 
of the scattering process (e.g. in DVCS it is related to the Bjorken variable). This leads to 
the so-called skewed distributions^, families of generalized parton distributions depending 
explicitly on the "skewedness" parameter, £ [^ |3|, f|, [J. In another approach, proposed 
by Radyushkin ||, [5| , one writes a spectral representation for the matrix element of the 
light-ray operator as an independent function of p ■ z and r ■ z in terms of a so-called 
double distribution. The skewed distribution for a given value of £ is then obtained as 
a particular one-dimensional reduction of this two-variable distribution. The advantage 
of this approach is that it allows one to make statements about the dependence of the 
skewed distribution on the skewedness parameter, £. 

The general structure of skewed and double distributions -- their symmetries, limit- 
ing cases, possible singularities, etc. - is a problem of great theoretical and practical 
importance. This problem has two aspects. The distributions depend, of course, on the 
behavior of the matrix elements of the light-ray operators as functions oip-z, r-z. This is 
a dynamical question, which one can address from the point of view of general invariance 

lr rhe term "skewed distribution" has been recommended as a common name for the "off-forward" 
distributions introduced by Ji @, |) and the "non- forward" distributions of Radyushkin |4|, ||, which 
differ in the definition of the parton momenta, see Rcf.p] for a detailed discussion. It encompasses 
also the "non-diagonal" distributions parametrizing matrix elements between hadron states of different 
quantum numbers, as have been introduced e.g. to describe DVCS with N-A transitions JllJ. 



principles, or by calculations using some dynamical model. However, the properties of the 
distributions are also determined by the particular way in which one writes the spectral 
representation for the matrix element. This concerns such things as e.g. the number of 
independent "twist-2 structures" one includes in the double distribution representation 
of the matrix element. A clear understanding of both aspects of this problem is necessary 
for building realistic models for skewed distributions. 

In this paper we investigate the structure of hadronic matrix elements of twist-2 oper- 
ators at a low normalization point, using general principles (symmetries, crossing etc.) as 
well as specific dynamical models, and consider the implications for spectral representa- 
tions in terms of skewed and double distributions. We show that the standard definition 
of the double distribution representation of Refs. ||, || is not always compatible with the 
basic structure of the matrix elements (to the very least, it implies severe singularities 
of the double distribution), and propose a complete representation which explicitly takes 
into account all twist-2 structures. The additional terms give rise to contributions to the 
skewed distribution which are non-zero only in the region — £/2 < X < £/2, and thus 
naturally lead to a "two-component" form of the skewed distribution, i.e., to essentially 
different functions in the two regions |X| < £/2 and \X\ > £/2. Such behavior was first 
observed in a model calculation of the flavor-singlet skewed distribution in the large-iV c 



limit in Ref.[12 



We find it useful to consider in addition to the nucleon matrix elements of twist- 
2 light-ray operators also the matrix elements between pion states. While hardly the 
target of choice for actual DVCS experiments, the pion is interesting from a theoretical 
point of view, for various reasons. First, it allows to avoid complications due to spin, and 
also its mass can be neglected. Second, the interactions of the pion with external fields are 
described completely by the chiral Lagrangian, which makes it possible to derive certain 
sum rules for the skewed distributions at a low normalization point from first principles. 
Finally, both the skewed and the double distribution in the pion at a low normalization 
point can be estimated in the large-iV c limit in the effective low-energy theory derived 



from the instanton vacuum of QCD |]13| . This is a fully field-theoretic description of the 
pion, which respects general properties such as crossing symmetry etc., and incorporates 
the consequences of the dynamical breaking of chiral symmetry. The same approach 
has been shown to give a realistic description of the quark/antiquark distributions in 
the nucleon (both usual [14 1 and skewed fL2| ) as well as the pion distribution amplitude 



JTR O]. In fact, the results of our calculation of the skewed and double distributions in 
the pion fully support our general conclusions concerning the need to modify the double 
distribution representation of Refs. |4|, [|] and the "two-component" structure of the skewed 
distribution.^ 

Another reason for our interest in the pion is the fact that the process related to DVCS 



2 In the case of the nucleon the calculation of non-forward matrix elements is complicated by the 
parametric restrictions imposed on the different components of the nucleons' momenta by the 1/N C - 
expansion. While it is possible to compute within the standard l/7V c -expansion the skewed distributions 
in the nucleon in the parametric range X, £ ~ 1/N C [12], it is difficult to get the double distribution in 
the nucleon in this approach; see Section |3|. 



off the pion by crossing, namely production of two pions in 7*7 collisions, can be mea- 
sured at low invariant masses [|17]. This process is described by a two-pion distribution 
amplitude, which, by crossing, is related to the skewed parton distribution in the pion 
fl8| , |I9| . In this paper we derive the explicit relation between the two functions, using 
dispersion relations to connect the regions of spacelike and timelike momentum transfers. 
In particular, this relation allows us to connect moments of the usual quark/ ant iquark 
distribution in the pion to characteristics of distribution amplitudes of two-pion reso- 
nances [[HJ. Given the contributions that measurements of 7*7 — ► 7r° |20| have made to 
our knowledge of the single-pion distribution amplitude, study of the process 7*7 — ► 2n 
could well be one of the cleanest ways to get information about the quark distributions 
in the pion — skewed as well as usual. 

The scale dependence of skewed and double distributions is described by generalized 
evolution equations, which combine features of both the DGLAP evolution for usual par- 
ton distributions and the Efremov-Radyushkin-Brodsky-Lepage evolution pTJ for meson 



distribution amplitudes. This problem has extensively been treated in the literature, see 
e.g. Refs. [§, |, |, |[ pi], |2|, |23|, |24j| . We shall not be concerned with this aspect here, but 
rather focus on the structure of the distributions at a low normalization point, how they 
are constrained by general principles (symmetries, crossing, etc.) and how they can be 
estimated in dynamical models taking into account non-perturbative effects such as the 
dynamical breaking of chiral symmetry, etc. 

We shall proceed as follows. In Section ^| we discuss the properties of non-forward 
hadronic matrix elements of twist-2 operators and their spectral representation from a 
general point of view. In Subsection |2.2|, using the pion as the simplest example, we show 



the importance of explicitly including all twist-2 structures in the double distribution 
representation, and discuss the implications for skewed distributions. The investigation 



is extended to nucleon matrix elements in Subsection |2.3| , with analogous conclusions. 
In Section || we perform a model calculation of the non-forward pion matrix elements 
and the corresponding skewed and double distributions at a low normalization point 
(/i ~ 600 MeV), using the effective low-energy theory based on the instanton vacuum. 
The results serve as an illustration for the general discussion in Section |2|. In Section [| 
we discuss the relation of the skewed distribution in the pion to the two-pion distribution 
amplitude. The crossing relation is derived in explicit form using moments. We use the 
crossing formula, together with the dispersion relation for the invariant-mass dependence 
of the two-pion distribution amplitude, to relate moments of the pion parton distribution 
to parameters of the distribution amplitudes of two-pion resonance wave functions. Our 
conclusions are summarized in Section [5| 

Appendix |A] gives a derivation of the generalized momentum sum rule for the skewed 
distributions in the pion. In Appendix [B] we consider "resonance exchange" contributions 
to the non-forward matrix elements in the pion. A general expression describing the 
contribution of the exchange of t-channel resonances of arbitrary spin is given. The 
results provide a simple dynamical explanation for the general properties of skewed and 
double distributions discussed in Section 0. 



2 Nonforward matrix elements and generalized par- 
ton distributions 

2.1 Skewed vs. double distributions: the pion 

To begin, we would like to discuss some general properties of non-forward hadronic matrix 
elements of QCD light-ray operators and their representation in terms of skewed and 
double distributions. We start with the simplest case, the pion, which already exhibits 
all features of interest to us here, and generalize to the nucleon in Subsection pT3[ 

Let us consider the non-forward matrix elements of twist-2 light-ray operators (nor- 
malized at some scale, /j,) between one-pion states. Due to isospin invariance, the matrix 
elements of the flavor-singlet and non-singlet quark operators are of the formP] 

<7r°(p - r/2)| H-z/2) I U b*/2, z/2] ^(z/2) \ix\p + r/2)> 

25 ab M I=0 (p-z, r-z; t) 

u r , W 

2ie abc M I=1 (p-z, r-z; t). 

Throughout the following the isoscalar and isovector parts of matrix elements in the pion 
will be understood to be defined as in Eq.([|); the isospin decomposition will often not be 
explicitly written. Here 1 and r 3 are flavor matrices; we consider the SU(2) flavor group. 
Furthermore, ip is the quark field, z^ a light-like distance (z 2 = 0), and z = 7 At z M . Finally, 
[—z/2, z/2] denotes the path-ordered exponential of the gauge field (phase factor) in the 
fundamental representation, 



[z/2, -z/2] = P exp 



1/2 

dt z^AJtz), 

1/2 



(2) 



which is required by gauge invariance; the path here is along the light-like direction, z. 
The matrix element of the corresponding twist-2 gluon operator is defined as 

^z»-r/2)| F, a (-z/2) l-z/2,z/2]F a u (z/2) \n b (p + r/2)) 

= 25 ab M G {p-z, r-z- t), (3) 

where F^ denotes the gluon field, and the phase factor is in the adjoint representation. 
In Eqs.(|H) and @ the dynamical information is contained in scalar functions, Ai I=0,1 



2 



and A4 G , which depend on the dimensionless invariants p-z and r-z, as well as of t = r 
From the mass shell conditions (p ± r/2) 2 = M 2 it follows that 

p-r = 0, p 2 = Ml- 1 -, (4) 



3 From Eq.(pJ) the matrix elements in charge eigenstates are obtained in the usual way: \ir°) = 
1 7r 3 ) , In^) = (In 1 ) ± i\-K 2 ))/^/2. Note that the neutral pion has no non-singlet matrix element due 
to C-invariance. 



so t is the only independent dimensionful invariant. In the physical region t < 0. In the 
following we consider the massless limit, M n — > 0. We note that G-parity (or, equivalently, 
time reversal invariance) requires that 

M I=0 (p-z, r-z) = M I=0 (-p-z,-r-z), 

M I=1 (p-z,r-z) = -M I=1 (-p-z,-r-z), (5) 

for same t; the symmetry of MP is the same as that of M I=0 . In fact, using in addition 
hermitean conjugation one obtains a stronger symmetry relating the functions with r-z — > 



—r-z and same p-z [25 



M I=0 ' l (p-z,r-z) = M I=0 '\p-z,-r-z); (6) 

this will be discussed in detail in Subsection fZ72 



Skewed distributions. In principle, the matrix elements Eqs.([|), (^) can be considered 
as functions of the invariants p-z and r-z as independent variables, defined in the physical 
region. However, in the amplitude for hard processes such as DVCS off the pion the 
matrix elements enter with some fixed ratio of r-z and p-z, 

r-z = £p-z, (7) 

which is dictated by the kinematics of the process; for instance, in DVCS £ is related to the 
Bjorken variable (— 1 < £/2 < 1) |2|, §, f|, ||. This suggests to define a "one-dimensional" 
spectral representation of the matrix elements in the form 



M 



7=0,1 



(p-z, r-z = £p-z; t) = 2p-z f dX e~ iXpz H I=0 '\X,£; t), 



where H I=0,l (X, £; t) are called the skewed quark distributions in the pion (the definition of 
the gluon distribution is analogous). The limits ±1 for the integral over the parameter X 
follow from rather general considerations ||, §, || . One can also give an explicit expression 
for H: Introducing a dimensionless light-like vector, n, and setting z = rn, one can invert 
Eq.(|) and obtains |, Q 



H{X,i;t) = i||V**-» 

x (tt(p - r/2)| i)(-rn/2) n [-tti/2, rn/2] ip(rn/2) \ir(p + r/2)) (9) 

[with isospin decomposition as in Eq. ([!])], and similarly for the gluon distribution. The 
symmetry property Eq.(^) requires H I=0 to be an odd function of X, H I=1 to be even. 
From the stronger symmetry Eq.(0) it follows that the skewed distribution is an even 



function of £ for any X; this was first noted in the case of the nucleon in Ref.p5[. 

The skewed distributions possess a simple partonic interpretation, the character of 
which depends on the relation of X to the skewedness, £; see Refs.0, |[ [26| for details. 
For X > £/2 and X < — £/2 the skewed quark distributions describe the amplitude for 



emission and reabsorption of a quark/antiquark in the infinite-momentum frame, and 
thus has properties analogous to the usual quark/antiquark distribution functions. For 
— £/2 < X < £/2, on the other hand, they have the character of distribution amplitudes 
for the creation of a quark/antiquark pair. One may thus expect the behavior of these 
functions to be quite different in the two regions. 

In particular, in the forward limit of the matrix element, r — > and £ — > 0, the skewed 
quark distributions reduce to the usual quark/antiquark distributions in the pion: 

H I=0 (X, £ = 0; t = 0) = ±[6(X)q s (X) - 6(-X)q s (-X)], 

H I=1 (X, £ = 0; t = 0) = 6(X)q v (X) + e(-X)q y (-X), (10) 

where q s (X),q v (X) correspond, respectively, to the singlet (quark plus antiquark) and 
valence (quark minus antiquark) distributions in a physical pion: 

q s {X) = [u + u] % ±{X) = [d + d]„±{X) = [u + u\^{X) = [d + d\„o(X), 
q v (X) = ±[u-u} 7T ±(X) = T [d-dU(X). (11) 

The moments of the skewed distribution, Eq.@, are given by non- forward matrix 
elements of local twist-2 spin-iV operators in the pion, which are parametrized by gener- 
alized form factors. On general grounds, the non-forward matrix elements of the spin-iV 
operators are irreducible rank-iV tensors constructed from the momenta p and r, so the 
moments of Eq.(^) are polynomials of degree at most N in £ || |26, 27]. In particular, 



the second moment of the isoscalar skewed distribution is related to the form factor of 
the QCD energy-momentum tensor ||. For the pion this form factor at t — can be 
computed from first principles using the chiral Lagrangian (see Appendix |S]), and one 
obtains a generalized momentum sum rule for the pion, 



dXX u '~" n<: 
-l 



H 1=U + H L 




(X,£;t = 0) = - 1 - V ■ (12 



The isovector skewed distribution in the pion is normalized to the pion electromagnetic 
form factor. For any t < 0: 

J^dXH'^(X^t) = F- m -(t). (13) 

Double distributions. Alternatively to the skewed distribution, Eq.(§), one can try to 
formulate a "two-dimensional" spectral representation of the matrix element Eq.fll]), as 
a function of r-z and p-z as independent variables. In the spirit of Refs.[(|, |5j we could 
write for the pion matrix elements a spectral representation in terms of a single function 
of two variables in the form 



M I=0 '\p-z, r-z; t) = 2p-z f dx e~ ixpz f % dy e - iyrz/2 F I=0 ' 1 {x,y; t), (14) 

J-i J-(l-\x\) 



where the functions F 1 °' 1 are called double distributions^. Here the range of the variables 
x, y is limited to j|, [5| 

-1 <x < 1, -(l-|x|) <y <l-|x|, 

see Fig.jl]. The property Eq.(|6]) implies that [p5| 



F^foy;*) = F^(x,-y;t). (15) 

The skewed distribution, Eq.(|]), is obtained as a one-dimensional "section" of this two- 
variable function, imposing a particular "skewedness" , £: 



f 1 dx I' W d^X-x-^F^x,^) = i^' 1 ^;*). (16) 

7-i 7-(i-|z|) 



This reduction process is illustrated in Fig.[l| In particular, in the forward limit the usual 
quark/ ant iquark distribution is recovered as 



r lXl dyF I =°(X,y;t = 0) = \ [9{X) q s (X) - $(-X) q s (-X)} (17) 

J — I1 — \X\) z 



and similarly for the isovector component, cf. Eq.(|lC|). 

The main reason for interest in a double distribution representation is the possibility to 
relate skewed distributions with different values of £. With certain assumptions about the 
behavior of the double distribution a number of statements about the skewed distributions 
follow immediately from the reduction formula, Eq. (|i6|) ; see Refs.[|], |5j for an extensive 
discussion. For instance, if the double distribution were continuous everywhere on its 
region of support (see Fig.|I|), the skewed distribution would be a continuous function of 
x and £. The double distribution is also convenient for model building, since any model 
of the double distribution, when inserted in the reduction formula, produces skewed dis- 
tributions satisfying the polynomiality condition for the moments |4], [5], [25|, ^6[. However, 



in order to be practically relevant, such applications require understanding of the general 
behavior of the double distributions, in particular, of their possible singularities. 

2.2 Trouble with double distributions 

When discussing properties of double distributions (such as their singularities) one should 
keep in mind that the behavior of these functions is determined by the behavior of the 
matrix element, Eq.(j^), as a function oip-z and r-z, as well as by the particular way 
in which one writes the spectral representation for it. This concerns, in particular, the 
number of independent "twist-2 structures" one takes into account in the decomposition 
of the matrix element. We shall argue now that it is not always adequate to represent the 
pion non-forward matrix element in the form of Eq. (0) , as a double spectral integral with 

4 We consider here the "modified" double distribution of Ref. S , which is appropriate for the symmetric 
choice of the momenta of the incoming and outgoing pion in Eq. (m) . 



a single prefactor, p-z. This form is incompatible with general features of the dependence 
of the isoscalar matrix element, Eq. ([!]), on p-z and r-z, and insisting on it one would incur 
severe singularities in the double distribution. 

In order to obtain information about the behavior of the function Ai(pz, r-z, t), Eq. ([!]), 
it is useful to consider instead of Eq.(|I|) the more general matrix element of the non-local 
vector operator [isospin components are defined in analogy to Eq.(|])] 

(-K(p-r/2)\$(-z/2)^[-z/2,z/2}iP(z/2)\ir(p + r/2)) = M„(p,r,z), (18) 

from which the matrix element Eq.(|l]) is obtained by contraction with the light-cone 
vector, z^. On general grounds, this matrix element can be parametrized as (for both 
I = and 1) 

M^(p,r,z) = 2p ll G + r M G|| + term oc z^, (19) 

where G and G\\ are generalized form factors depending on p-z,r-z and t. The terms 
proportional to z^ vanish upon contraction with z^ and does not contribute to the twist-2 
part matrix element Eq. ([]]). The term proportional to r M , however, does contribute to 
Eq.([ip. In fact, it is the presence of this structure which causes trouble in the double- 
distribution representation of the isoscalar matrix element, Eq . (|14D . 

In the limit £ — > the operator in Eq.(^) reduces to the local vector current, which 
is conserved. This implies that G\\(z = 0) = for all t, i.e., the matrix element is "trans- 
verse" (oc p^). However, the non-local operator with z ^ is generally not conserved, so 
there is no reason for G\\ to be zero for z ^ 0. Actually, current conservation is only a 
sufficient condition for G\\ to be zero, not a necessary one. For z = one obtains Gu = 
already from time reversal invariance and the hermiticity of the local current operator. 
Applying the same symmetry transformations to the non-local operator, one finds that 
for the isoscalar matrix element 

G I=0 {p-z, r-z, t) = G I=0 {p-z, -r-z, t), 

Gl=°(p-z, r-z, t) = -G^ip-z, -r-z, t). (20) 

In the local case r-z = 0, and Gl =0 would be zero identically in t. However, in the general 
case, 2^0, there is again no reason for Gjj =0 to be zero. 

The presence of a "longitudinal" (oc r M ) part of the vector matrix element, Eq. (|T9"D , 
means that the matrix element M. I=0 , obtained by contracting Eq . ([I9|) with z M , contains 
in addition to the p-z-teim. a piece with prefactor r-z, 

M I=0 {p-z,r-z,t) = p-zG I=0 {p-z,r-z,t) + r-zG\^{p-z,r-z,t). (21) 

In particular, since generally Gl =0 ^ 0, M does not vanish in the limit p-z — > and 

r-z 7^ 0: 

M I=0 (p-z^0,r-z^0,t) ^ 0. (22) 

10 



One can easily see that this implies that the matrix element Ai I=0 cannot be represented in 
the form Eq.flHD with a non-singular double distribution F I=0 (x,y). Suppose F I=0 (x,y) 
were non-singular in its range of support. It would then define the matrix element on the 
L.H.S. as an analytic function of the variables p-z and r-z, which could be continued to a 
point (in the unphysical region) where p-z = 0, but r-z/0. At this point the R.H.S. of 
Eq.(|14"D vanishes because of the prefactor p-z, but not the matrix element on the L.H.S. , 
cf. Eg. fl2"2"|) . Clearly, this implies that F I=0 (x,y) must be singular, in one way or another. 
What would be the character of these singularities? Trying to absorb a pz-independent 
piece in the integral Eq. fllil) would amount to finding an integral representation of I /p-z 
in the form 

— = f 1 dxe~ lxpz f(x), (23) 

p-z J-i 

with f(x) some generalized function. Assuming that the integral on the R.H.S. can be 
continued to p-z — > 0, one concludes that no Mellin moments of the function f(x) exist. In 
particular, this means that the singularity in f(x) cannot be of delta-function type. (We 
shall return to this point below.) Thus, we conclude that, although the two contributions 
to M. I=0 in Eq.fl2"T|) are not structurally distinct, it is not possible to include the r-z-term 
in a double distribution representation of the form Eq.(|14D staying within the usual class 
of generalized functions. 

It is important to mention that the difficulties noted here do not concern the repre- 
sentation of the matrix element through a skewed distribution, Eq.(||]). In this case one 
writes the representation of the matrix element under the condition r-z = £,p-z, with £ 
fixed. In particular, repeating the above argument and taking in Eq.(||) the limit p-z — > 
we would now also have r-z — > 0, so that both L.H.S. and R.H.S. of Eq.(||) would vanish, 
and there is no need for H I=0 to be singular. 

Since there are many advantages in a double-distribution representation of the non- 
forward matrix elements, it is worthwhile to think how Eq.(|14D could be modified to allow 
for a double spectral representation in terms of standard generalized function. The origin 
of the problems with the form Eq.fll^D is that the matrix element does not go to zero in 



the limit p-z — > 0, Eq . (|22j) . One possibility would be to simply omit the prefactor p-z 
in Eq. (|I~4"D ; however, this would result in a function F I=0 (x,y) which does not reduce to 
the usual parton distribution in the forward limit, r — > 0, cf. Eq . (|T7|) , and would not be 
useful for model building. Alternatively, one could add to Eq.(|TJ]) a term not vanishing 
in the limit p-z — > 0. Minimally, this could be a term depending only on r-z, which can 
be represented by a one-dimensional spectral integral: 



M I=0 {p-z, r-z;t) = 2p-z f dx e~ ixpz f dy e~ iyrz/2 F I=0 {x,y;t) 

J-l J-ll-\x\) 



»l-\x\ 

'-(1-1*1) 

+ r-z f 1 dye- iyrz/2 D{y;t). (24) 



The functions F I ~°(x, y; t) and D(y; t) are uniquely defined if we understand the first term 
to be a representation of M. I=0 {p-z, r-z; t) — .M /=0 (0, r-z; t), the second of .M /=0 (0, r-z; t), 
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i.e., as a "subtraction term". The explicit factor r-z in front of the second term is natural 
since J\A I=0 (0,r- z — ► 0) = 0. The support of D is limited to — 1 < y < 1, i.e., this 
function has the character of a distribution amplitude. Time reversal and hermiticity, 
Eg. fl20|) , require that (same t) 

F /=0 (x,y) = F I=0 (x,-y), D(y) = -D(-y), (25) 

i.e., the behavior with respect to y — > — y of the new function D(y) is opposite to that of 
the usual double distribution, Eq.(|l5|). 

The skewed distribution which follows from the new representation Eq.([24j) is now the 
sum of two contributions: 



H I=0 {X,&t) = I' dx /* '*' dy6(X-x-y^2)F I=0 (x,y;t) 

J-l J-(l-\x\) 



-1 fl-\x 

-(1-1*1) 

+ sign^feiV (26) 

Note that both contributions are even functions of £, in accordance with the general 
symmetry of the skewed distribution following from Eq.(^). The first piece follows the 
usual reduction formula, Eq.(^), and is generally non-zero in the entire range — 1 < X < 
1. The second piece is obtained by substituting in Eq . fl2"4|) r-z = £,p-z and changing 
the integration variable to £y/2. Since the support of D is limited to — 1 < y < 1 this 
contribution to the skewed distribution is present only for — £/2 < X < £/2. Thus, the 
need to include the "subtraction term" in the double distribution representation naturally 
leads to a skewed distribution with essentially different behavior in the regions \X\ < £/2 
and in \X\ > f/2. 

So far we have explored the consequences of the presence of "longitudinal" terms in 
the matrix element Eq.(|l9|), or of property Eq. (p2|) , from a general point of view, arguing 
that there is no reason for such contribution to be zero. In Section [5] we perform a model 
calculation of the pion matrix elements at a low normalization point in the effective chiral 
theory based on the instanton vacuum, wich shows that such terms in the matrix element 
do indeed appear, and lead to the "two-component" form of the skewed distribution 
described above. 

The region p-z — > and r-z ^ implied in the limit in Eq. (|22| ) corresponds to values of 
|£| > 2, which are not physically accessible in DVCS. However, using crossing invariance 
one can relate the function Ai(p-z,r-z) in the unphysical region, Eq.fl2"2"|), to the matrix 
element for two-pion production by a light-ray operator in the physical region (see Section 
[|). The latter can be measured (e.g. in 7*7 — > nn and 7*iV — > ttttN reactions) and is 
generally non-zero, providing additional evidence for the presence of r-z-terms and the 
property Eq. ([Z2"D . 

A simple dynamical explanation for the origin of r-z-terms in the isoscalar pion matrix 
element can be found by considering "resonance exchange" contributions to the matrix 
element, in the spirit of the vector dominance model for the pion electromagnetic form 



factor |3q| . By this we mean "factorized" contributions to the matrix element in which 
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the pion and the light-ray operator communicate by t-channel exchange of a resonance 
characterized by a twist-2 distribution amplitude, as are shown schematically in Fig.Q. 
In Appendix [B] we derive a general formula describing the contribution resulting from the 
exchange of a resonance of arbitrary spin to the pion matrix element. In particular, we 
show there that the property Eq. ( [2"2"D of the isoscalar matrix element is naturally obtained 
from exchange of even-spin isoscalar resonances. Far from being a complete dynamical 
description of the non- forward matrix element, this phenomeno logical model helps to 
develop an intuitive understanding why the structures described above appear. 

The double distribution F I=0 (x,y) in Eq . (f2"4|) is a generalized function which may 
contain delta function type singularities. Such terms in the double distribution [in the 



restricted ansatz Eq. fll4|) 1 were studied by Radyushkin in connection with resonance ex- 
change contributions to the non-forward nucleon matrix elements j5], |2B[. We already 
argued above that the pure r-z-terms in the isoscalar pion matrix element, cf. Eq.(|22|), 
cannot be described by delta function type contributions to F I=0 (x,y). It is interest- 
ing to verify this at the level of the reduction formula, Eq. (|26|) . Could delta function 
contributions to F I=0 (x,y) mock up the structure of the D-term in H I=0 7 A term in 
F I=0 (x,y) of the form 5(x)(j)(y) would give a contribution to the skewed distribution 
H I=0 oc <f)(2X/£) /|£|, which is ruled out because H I=0 must at the same time be odd in 
X and even in £. One thus has to turn to derivatives of 5(x). A term S'(x)x(y) would give 
a contribution H I=0 oc sign(£) x'(2X/£) /£ 2 , which can be non-zero but cannot describe 
the contribution generated by D(y) in Eq.(f26|) (consider for example the forward limit). 
This argument can easily be extended to any derivative of 5(x). Thus, we conclude that 
the term generated by D(y) in Eq.fl2~3]) represents a genuine separate structure which 
cannot be obtained from delta function contributions to F I=0 (x, y)^\ 

We remark that in the resonance exchange model of Appendix [FJ exchange of even-spin 
resonances generally contributes to both D(y) and to delta function terms in F I=0 (x,y), 
cf. Eq. flB.4Q . Spin-0 ("sigma meson") exchange is special in that it contributes only to 

D(y). 

In the amplitude for hard exclusive processes such as DVCS the skewed distribution 
is convoluted with a hard scattering kernel which is singular at X = ±£/2 0, S, |j, 0]. 
For this integral to exist (i.e., for factorization to hold) it is important that the skewed 
distribution be continuous in X at these points. Assuming that the first term on the 



R.H.S. of Eq.(p6|) is continuous, this would be satisfied if 

D(y-t) - (y->±l). (27) 

We shall see below that it is indeed reasonable to expect that D(y; t) satisfies this property, 
reminiscent of a meson distribution amplitude. Model calculations of the matrix elements 
in the effective chiral theory bases on the instanton vacuum and in a "resonance exchange" 
model give rise to functions satisfying Eq . (|27|) . 

The representation Eq.(|2~3]) is the minimal modification of Eq.([TJ|) consistent with 
Eq . (|22"D . For some purposes (e.g. crossing symmetry) it could be convenient to have a 

5 It is amusing to note that, rather than a derivative of 8(x) this term represents, in a sense, an 
"integral" of 6(x), cf. Eq.@. 
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representation which is symmetric with respect to p-z and r-z. One could write: 
M I=0 (p-z, r-z) = 2p-z f dxe~ ixpz [ ~* dy e~ iyrz/2 F I=0 (x,y) 

J-l J-(l-\x\) 



(1-1*1) 

+ r-z I dxe~ ixpz f X dy e~ iyrz/2 Ft =0 (x,y). (28) 
1 ■'-(1-1*1) 



The two functions F 1 ° and F* ° would be uniquely determined if we defined them as the 

r ll 



spectral representation of the generalized form factors G 1 °{p-z, r-z, t) and Gl °(p-z, r-z, t) 



of the vector operator, Eq.([L9]). 6 Again, Eq.([20"D requires that 

F'=°(x,y) = F^^-y), F{=°(x,y) = -F{=°(x,-y). (29) 

In terms of these new functions the skewed distribution would now be given by the re- 
duction formula Eq.(|T6|) with 

F(x,y) -> F(x,y) + ^F\\(x,y). (30) 

Note that both F and F» are generalized functions which may contain delta-function 
singularities. 

Finally, let us note that for the isovector pion matrix element, Ai I=1 , the original form 
of the double distribution representation, Eq.(|TJj), works fine, and no "subtraction terms" 
of the kind in Eq.(|24D are required. This is because M. I=1 is odd in p-z for any r-z and t 
[as follows from combining Eq.([5[) and Eq.(|6])], and thus 

M I=1 (p-z^0,r-z^0,t) -> 0. (31) 

In the resonance exchange model this is again easily understood; it follows from the fact 
that isovector two-pion resonances have odd spin, see Eq.( |B.4| ) in Appendix [l|. 

2.3 The nucleon 

We now turn to non-forward matrix elements in the nucleon. By a simple extension of the 
arguments offered in the previous subsection for the pion, we show that also in the case 
of the nucleon the double distribution representation of Refs. 0, |5|, ^B| should be modified 
to take into account all possible twist-2 structures. 

The object of interest now is the nucleon matrix element of the twist-2 light ray 
operator of Eq. ([]]). Again we distinguish the isoscalar and isovector matrix elements: 

(N(p - r/2)X, A' | 1>{-z/2) \ 1 z [-z/2, z/2] ^(z/2) \N(p + r/2), T 3 , A) 



6 Strictly speaking, we have no general proof that a double spectral representation for the form factors 
G and G» exists. At least in our model calculations in Section 13 and Appendix pi we shall encounter 
only contributions to the matrix element which can be represented by Eq.pH) with F, Fu having at most 
delta function singularities. 

14 



28 Tan M I=0 (\',\;p-z } r-z; t) 
2(T% iT ,M I=1 (\',\;p-z, r-z; t), 



(32) 



where T^,T^ denote the isospin projection (Ts,T^ = ±1/2 for proton/neutron). The only 
difference to the pion is that now the functions Ai I=0,1 depend also on the helicities 
of the incoming (A) and outgoing (A') nucleon. In analogy to the pion, Eq. QT9|) , let us 
consider also the matrix element of the more basic light-ray operator with 7^ [the isospin 
decomposition is analogous to Eq. (|32]) and not written], 

(N(p-r/2),X'\i,(-z/2) ltl [~z/2,z/2]^z/2) \N(p + r/2),X) = M,(X\ X;p,r, z), 

(33) 

which on general grounds can be parametrized as (for both 1 = and 1) 



M L 



U' 



7^1 + 



2M 



°uvr u G 2 + r u G\\ + ... 



N 



U 



(34) 



where U' = U(p — r/2, A'), U = U(p + r/2, A) are the nucleon spinors. We have not 
written explicitly terms which vanish upon contraction with z^ and do not contribute 
to the twist-2 part (such as z^, a^ v z u ). Here Gi,G 2 ,G\\ are generalized form factors 
depending on p-z,r-z and t. In the limit of a local operator, z — > 0, Gi and G 2 reduce 
to the usual Dirac form factors for the vector current, and G\\ — > because of current 
conservation. However, as in the case of the pion, for z ^ the term oc r^ is generally 
non-zero, G\\ 7^ 0. 

Contracting Eq. (|34]) with z^ we obtain the twist-2 matrix element, Eq. (|32]) (for both 
/ = and 1): 



M(X',X;p-z, r-z;t) 



U' 



2G ± + ^ir a ^ rVG i + r ' zG \\ 

Zivijq 



u. 



(35) 



Actually, here the three structures are not independent. Using the well-known identity 

U'Upp = 2M N U' ltl U + iU'a^Ur", (36) 

we could rewrite the third term as a linear combination of the first two. In this way we 
would arrive at (I = 0, 1): 



where 



M{X',X;p-z, r-z;t) 



~ V ' z 

G\ = Cti t" 2 ^11) 

p-z 



U' 



zG x 



2M N ^ 



o,,„z»r v G, 



U, 



Go 



G 2 + 2—G\\. 
p-z 



(37) 



(3* 
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A decomposition of the form Eq.(^) was assumed in Refs.[|], |j, where a double distri- 
bution representation of the matrix element was proposed in the form 

M I=0 >\\',\;p-z, r-z-t) 
= U'zU [ dxe~ ixpz [ X dy e~ iyre/2 F I=0 >Hx,y\t) 

J-l y_(i_|a;|) 

+ -^z»r u U'a, LU U [ dxe~ lxpz f ~* dy e~ iyrz/2 K I=0 Hx,y,t). (39) 
2M N ^ J -i J-{i-[x\) 

We see that in the isoscalar case this ansatz suffers from the same problem as the simple 
ansatz for the double distribution in the pion, Eg. fll4|) . Since in general Gj| =0 7^ 0, the 

functions G{ =0 ,G2 =0 in the reduced decomposition, Eq . (|38|) , have singularities of the type 
1/p-z, which cannot be represented by spectral integrals with usual generalized functions, 
as described in the previous subsection. Thus, the conclusion is the same as for the 
pion: Although the twist-2 contributions from the "longitudinal" part (oc r M ) of the 
vector matrix element are not structurally distinct from those from the "transverse" part 
(oc p/j), one cannot obtain them from a double distribution representation whose form is 
modeled on the "transverse" part. 

Again we stress that there is no problem with a representation of the matrix element 
Eq.(|3~?D in terms of skewed distributions (see Refs.0, |3|, [|, Q for their definition in the 
nucleon). In this case the factors r-z/p-z incurred in eliminating the r-z-term, Eq. fl3"8f ), 
are replaced by the skewedness, £, which is a fixed external parameter. 

In analogy to the pion, Eq.([24|), we suggest to modify the spectral representation 
Eq. ( |3"9"D by explicitly including the r • z-terms. A minimal variant would be to add a 
term depending only on r-z, in which one could absorb the p ■ ^-independent part of the 
"longitudinal" term of Eq. fl2ll) , r-zGl =0 (p-z = 0,r-z): 

M I=0 (\',\;p-z, r-z-t) 
= U'zU f dxe- ixpz ( l ~ X dye- iyrz/2 F I=0 (x,y;t) 

J-l J-(l-\x\) 

+ ^— z » r v\j> 0ilv U f 1 dxe~ ixpz j 1 ' X dye- tyrz/2 K I=0 (x,y;t) 
2M N ^ J-i J-(i-\x\) K J 

+ U'Ur-z C dy e- iyrz/2 D(y;t). (40) 

Alternatively, one could directly work with the spectral representation of the form factors 
Gi,G 2 and G\\, as in the representation Eq. (E3) for the pion. 

The skewed distribution in the nucleon resulting from the full double distribution rep- 
resentation, Eq. (|40|) , again has a "two-component" form, since the term with D(y) gives 
rise to a contribution non-zero only in the region — £/2 < X < £/2. (The correspond- 
ing reduction formulas can be obtained by a trivial modification of the ones written in 
Refs.[[|, pi p5[.) This explains the behavior of the isoscalar skewed distribution, H I=0 , 
which was encountered in a model calculation in the large- iV c limit [I2~|. 
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3 Distributions in the pion from effective chiral dy- 
namics 

For quantitative estimates of the non- forward matrix elements Eq.(|J) and the skewed 
and double distributions one has to turn to model calculations. Here we compute these 
quantities at a low normalization point in the low-energy effective field theory based on the 
instanton model of the QCD vacuum. This effective theory incorporates the dynamical 
breaking of chiral symmetry, and provides a realistic description of hadronic properties 
of the pion and nucleon |13, [2£|. Its content can be summarized in an effective action 
describing the interaction of a pion field with massive "constituent" quarks, in a way 
which is dictated by chiral invariance: 

S cS = fd 4 X ij(x) ltd - MF(d 2 ) e^r-ff-W/F, F ( d 2j\ ^ x y (41) 



Here, 7r a is the pion field, and F n = 93 MeV is the weak pion decay constant. The dynam- 
ical quark mass generated in the spontaneous breaking of chiral symmetry is momentum 
dependent; the form factors F(d 2 ) are related to the instanton zero modes |13||. They cut 
loop integrals at momenta of order of the inverse average instanton size, p _1 ~ 600 MeV. 

The effective theory Eq.([|l]) has been derived from the instanton model of the QCD 
vacuum. This allows for an unambiguous identification of the twist -2 QCD operators with 
operators in the effective theory. It is understood that the QCD operators are normalized 
at a scale of the order p = p _1 rs 600 MeV. The general framework for computing 
parton distributions and meson wave functions at a low normalization point within this 
approach has been developed in Refs. |14], |T5|1 . An essential point is that the value of the 
dynamical quark mass, M, is parametrically small compared to the UV cutoff, p _1 ; their 
ratio is proportional to the packing fraction of the instanton medium (Mp) 2 ~ (p/R) 4 . 
Qualitatively speaking this means that in leading order in this parameter one is dealing 
with structureless constituent quarks; in particular, the gluon distribution appears only 
at order (Mp) 2 . At a technical level, working in leading order in Mp means retaining 
only the ultraviolet divergent part of the quark loop integrals computed with Eq . (|4lD , 
absorbing the ultraviolet divergence in the pion decay constant, F n . 

The non-linear form of the coupling of the pion to the quarks in Eq.(f|l]) is required 
by chiral invariance. Expanding the exponential in powers of the pion field we obtain 

The effective theory contains a Yukawa-type quark-pion vertex as well as a two-pion 
quark vertex. Consequently, there are in general two contributions to the matrix element 
of a twist-2 quark operator between pion states, corresponding to the diagrams (a) and 
(b) of Fig.[3[ The diagram (a) of Fig.|3] contributes only to the flavor-singlet matrix ele- 
ment, while (b) contributes both in the flavor-singlet and non-singlet case. The Feynman 
integrals can straightforwardly be computed introducing light-cone coordinates with re- 
spect to the vector n ~ z; see Refs.[|T3], |T^, |TB| for details. The integral over transverse 
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momenta contains a logarithmic divergence which is cut by the form factors, F(d 2 ). More 
simply, one may keep only the logarithmically divergent part of the diagram and absorb 
the logarithmic divergence in the pion decay constant. It was shown in Refs.fO, Pj] that 
this is a legitimate approximation except in the vicinity X = ±£/2. (We shall include the 
form factors in the calculation later.) 

In this approximation the contributions of diagrams (a) and (b) of Fig.|3] to the matrix 
elements, Eq.([TJ), can be computed analytically. For the isoscalar part we obtain (for 
simplicity we take t — > 0): 



M I=0 (p-z,r-z) {a) = 2% 
M I=0 (p-z,r-z) (b) = 2i 



r-z 2 r-z 

cos 1 sin 

2 r-z 2 



2 . r-z 
cos p-z sin 

r-z 2 



and the total result is 



M I ~°(p-z,r-z) 



2i 



r-z 



cosp-z — cos ■ 



(43) 
(44) 

(45) 



Here contribution (a) depends only on r-z; due to the contact nature of the two-pion-quark 
vertex, Eq. (pE2"D , the average momentum p does not enter in the quark loop, see Fig.[3]. This 
contribution vanishes in the forward limit r — > 0. Note that both contributions to the 
isoscalar matrix element, as well as their total, behave as described in Section ||]: they do 
not go to zero in the limit p-z — > 0, r-z ^ 0, and thus cannot be represented by a double 
distribution in the form Eq. fll4]) . Within the proposed new representation, Eq. (|24D , which 
allows for a ^-^-independent part, this model result would correspond to 



F I=0 (x,y) 

D(y) 



- [0(-l < x < 0) - 6(0 < x < 1)] 6{y), 
i[-0(-l<y<O) + 9(0 <y<l)}, 



(46) 
(47) 



where 9(a < X < b) is unity if a < X < b and else zero. That F I=0 here is proportional 
to a delta function in y should be seen as an artifact of keeping only the logarithmically 
divergent piece; this would change when retaining finite terms at t ^ 0. 
The result for the isovector matrix element, Eq. ([!]), at t — > is 



M 1 1 (p-z,r-z) = M 1 1 (p-z,r-z) 



(b) 



4 sin p-z. 



(48) 



This matrix element vanishes in the limit p-z —>■ 0, hence there is no problem with 
representing it by a double distribution in the form Eq. (|H|) : 



?i=i. 



x,y) 



0(-l <x< l)S(y). 



(49) 



The corresponding skewed distributions may be computed either using the results for 
the matrix elements, Eqs.(H3T) and (|4"4"[), and the definition Eq.(|8l), or directly by computing 
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the R.H.S. of Eq.(||); both ways lead to identical results. For the isoscalar part we find 

# /=0 (X,0 (a) = -~ <H-t;<X <^\^, (50) 



H I=0 (X,0 



(b) 




-6 [-KX <- 



H\^<X<1 




the total result is 



H 1=0 {X,£) 



•1< X < 



<X < 1 



(51) 



(52) 



The functions are shown in Fig.|j. One sees that the contribution from diagram (a), 
Eq . (|50"D , is non-zero only in the region — £/2 < X < £/2. It is absent in the forward 
limit, r — » 0. Note that this contribution to the skewed distribution is discontinuous 
in X at ±£/2; this behavior will be modified when taking into account the momentum 
dependence of the dynamical quark mass, see below. The contribution from diagram (b) 
is continuous at X = ±£/2; this part reduces in the forward limit, £ — > 0, to the singlet 
quark distribution in the pion, Eq.(ID), which in this approximation (keeping only the 
logarithmic divergence, neglecting the form factors) would simply be given by 



Qsi 



sing 



X 



6(0 <X < 1). 



(53) 



The result Eq.([52"D is consistent with the generalized momentum sum rule, Eq. fll2"|) . 
In our approach based on the instanton vacuum the gluon distribution is parametrically 
small, oc (Mp) 2 ~ (p/R) 4 , so the skewed quark distribution, Eq.(^), should saturate the 
sum rule at the low normalization point. Integrating Eq.([52|) we observe that, indeed, 



fdXXH'^iX^) = Ul 



(54) 



Eqs. (|50|) , fl5"l~|) and ([52] ) represent the result for the skewed distribution obtained with- 
out taking into account the momentum dependence of the dynamical quark mass. The 
discontinuity at X = ±£/2 in the contribution (a) to H I=0 obtained in this approxima- 
tion would violate the factorization of the DVCS amplitude, since the hard scattering 
kernel contains poles at X = ±£/2. However, as was shown in Refs.[T^], the momentum 
dependence of the dynamical quark mass can not be neglected for values of X near ±£/2, 
since in this case the integral over transverse momenta is cut by the form factors already 
at momenta of order M ^C p _1 . The same mechanism makes the pion distribution ampli- 
tude vanish at the end points [[HJ. In Fig.|] we show the two contributions (a) and (b) to 
H I=0 obtained when taking into account the form factors F(d 2 ) (we use the simple ana- 
lytic approximation of Eq.(24) of Ref.||15|). As expected, contribution (a) now vanishes 
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at X = ±£/2, while the modification of contribution (b), which was continuous already 
without form factors, is only quantitative. 

Finally, the result for the isovector skewed distribution at t — > is 

H I= \X,0 = tf /=1 (X,0 (b) = 9(-KX<l), (55) 

which is an even function of X, in agreement with C-invariance. [In the forward limit 
this corresponds to a valence quark distribution in the pion q v (X) = 9 (0 < X < 1), cf. 
Eq . fllPf ) , so comparing with Eq.(|53D we see that in this approximation (no form factors) 
the "sea" quark distribution in the pion is zero.] As in the case of contribution (b) to the 
isoscalar distribution, inclusion of the form factors F(d 2 ) does not change the result for 
the isovector skewed distribution in an essential way, except for forcing the distribution 
to vanish at X — ±1.Q 

Some comments are in order concerning the calculation of skewed and double distribu- 
tions in the nucleon. In the large-iV c limit the nucleon in the effective low-energy theory 
is characterized by a classical pion field ("soliton") |28|| . Quantization of the translational 
and rotational zero modes in the framework of the l/iV c -expansion gives rise to nucleon 
states with definite momentum and spin/isospin quantum numbers. When applying this 
approach to the computation of non-forward matrix elements of the type Eq.(|l]), the stan- 
dard l/iV c -expansion implies that different components of the average momentum, p, and 
momentum transfer, r, are of different order in N c [the nucleon mass is 0(N C ), while the 
momentum transfer in the Breit frame is 0(N®)}. While it is possible to compute the 
skewed distribution, Eq.(|^), in the parametric range X, £ ~ 1/N C JT2J, it is not possible 



to uniformly obtain the matrix element Eq.([|) in the whole kinematical range necessary 
to restore the double distribution. In contrast, in the case of the pion all components of 
p and r are 0(N®), making it possible to treat p and r on the same footing. 

4 Crossing and the two— pion distribution amplitude 

4.1 Two— pion distribution amplitude 

An interesting feature of the pion is the fact that the quantity related to the skewed 
parton distribution by crossing, namely the two-pion distribution amplitude (DA), can be 
measured in two-pion production at low invariant masses. These 27rDA's were introduced 



recently in the context of the QCD description of the process 7*7 — ► 2n [|17]. We now 
establish explicitly their relation to the skewed quark distributions in the pion. This will 
allow us to express the ^-dependence of the lowest moments of the skewed distribution in 
terms of form factors in the timelike region. 

The two-pion DA's are defined, in analogy to the skewed parton distribution, as the 



The valence quark distribution in the pion has also been studied in the instanton vacuum in a 
somewhat different approach by Dorokhov and Tomio p9j. 
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matrix elements of the twist-2 operators between the vacuum and a two-pion state: 

ont (n a ( Pl )n\p 2 )\ $(-z/2) | c 1 z [-z/2, z/2] t[>(z/2) |0> 

,1 f 25 ab $ /=0 (u, C, W 2 ) 

= P-z du e i{u - 1/2)p - z { (56) 

Jo [ 2ie abc <5> I=1 {u,(,W 2 ). 

The outgoing pions have momenta pi,P2, and P = p\ + pi is the total momentum of the 
final state. The generalized DA's, Eq. (pB|), depend on the following kinematical variables: 
the quark momentum fraction with respect to the total momentum of the two-pion state, 
z; the variable ( = (z-pi)/(z-P) characterizing the distribution of longitudinal momentum 
between the two pions, and the invariant mass of the two-pion system, W 2 = P 2 . Also, 
an explicit representation of the DA can be written in analogy to Eq.(|J) [the isospin 
decomposition is analogous to Eq.([56|)1 

Xout(7r(pi)7r(p 2 )| ^(-rn/2) n [-Tn/2,rn/2] ^(rn/2) |0), (57) 

where n oc z is a dimensionless light-like vector. 

From C-parity one derives the following symmetry properties (we do not write the 
argument W 2 ): 

^=°(uX) = -$ /=0 (i-«,C) = $ /=0 (M-0, 

$ /=1 (n,C) = $ /=1 (l- M ,C) = -$ /=1 (n,l-C). (58) 

The first moment of the isovector (1 = 1) two-pion DA is the pion e.m. form factor in 
the time-like region, 

f 1 du<5> I=1 (u,(,W 2 ) = (2C- I) F°- m iW 2 ) , (59) 

Jo 

and thus scale-independent (F^- m -(0) = 1). For the isoscalar (/ = 0) part, however, we 
have the normalization condition 



f l du{2u-l)^=\uX,W 2 ) = -2M 2 W C(l-C)if MT (^ 2 ), (60) 

Jo 

where M 2 is the momentum fraction carried by quarks in the pion at the given scale, 
and F^ UT (W 2 ) is the form factor of the quark part of the energy momentum tensor, 
normalized to _F™ T (0) = 1. In Ref. |19[] this form factor was estimated in the instanton 
model of the QCD vacuum at low two-pion invariant mass: 

?v w 2 

F EMT (W 2 ) = 1 + c + . . . . 
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It is useful to expand the two-pion DA simultaneously in eigenfunctions of the ERBL 
evolution equation |H| [Gegenbauer polynomials Cf/ 2 (2u — 1)] and in partial waves of the 

1 /9 

produced two-pion system [Legendre polynomials, or Gegenbauer polynomials C l (2£ — 
1)]. Generically this decomposition is of the form [lTJ|: 



oo n+1 



HuX,W 2 ) = 6u(l-u)Y,Y, B m(W 2 )C s J 2 (2u-l)Cr(2C-l) 



(61) 



n=0 1=0 



where n runs over even (odd) and / over odd (even) integers for the isovector (isoscalar) 
DA, cf. Eq. Q5£|). The normalization condition Eq. fl5TJ| ) requires that B^iW 2 ) = 
F^ m '{W 2 ). Note that the asymptotic form of the isovector two-pion DA is given by |1 



$ 



i=i 

asymp 



u,(,W> 



6u(l-u)(2(-l)F°- m iW 2 



(62) 



In Ref. [|19| certain soft-pion theorems for the two-pion DA were proven, which apply 
in the regions ( — > or ( — > 1 and W 2 — > 0, where one of the produced pions becomes 
soft. In the isovector case (1 = 1) they relate the two-pion DA to the DA of one pion, 

4> n {u): 



<S> I=1 (u,( = 1,W 2 = 0) = -$ /=1 ( M ,C = 0,W 2 = 0) 
while in the the isoscalar case (I = 0) one obtains 

$ /=0 (u, C = 1, W 2 = 0) = <5> I=0 (u,( = 0,W 2 = 0) 



u 



(63) 



(64) 



The theorem Eq. (|63|) allows to relate the expansion coefficients of the isovector two-pion 
DA, Eq. (§T|), at W 2 = 0, with those of the pion DA, 



4>n(u) = 6u(l — U) 

the relation takes the form 
4.2 Crossing relation 



1 + E aWC^(2u-i; 



n+1 

Y: BifiW 2 = 0). 
1=1 



(65) 



(66) 



By crossing, the matrix element defining the two-pion DA, Eq.(|57|), is related to the one 
appearing in the definition of the skewed distribution, Eq.(^|). This allows one to express 
the moments of the skewed parton distribution in terms of the expansion coefficients of 
the two-pion DA, Eq. (pT]) . This relation takes the form: 



H N-l n+1 /f\ N 

/ dXX N ^H\X^t) = E Y,BUt) I c} 12 

J ^ 1 n=0 1=0 \ Z J 



A, 



x du 6u(l - u) {2u - if- 1 Cl /2 {2u - 1). (67) 
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One immediately notes that the R.H.S. is a polynomial of degree (at most) N in £, i.e., 
the polynomiality condition for the skewed distribution (see Subsection ^J] and Ref . ]27j ) 
is satisfied. Since only the B n i with n odd (even) are non-zero in the isoscalar (isovector) 
case, the skewed distribution is an odd (even) function in X, in agreement with G-parity, 
Eq.(|5[). Also, note that due to the restrictions in the values of I, cf. Eq.flBIp, H is an even 
function of £ for both 1 = and 1, as it should be. 

To prove the relation Eg. QS7D , we consider the expression for the N— th moment of the 
skewed distribution as a non-forward matrix element of a local spin-iV, twist-2 operator, 

f dXX N ~ x H(X,£,t) 

= (2p-n)' N (?r(p + r/2)|-0n-7 (n-vf'Vl^-^)). (68) 



The iV-th moments of the two-pion DA, Eg. (|57j) , is given by the vacuum to two-pion 
matrix element of the same local operator. Substituting the double expansion, Eg. flB"T|) , 
we obtain 



f 1 du(2u-l) N - 1 ^(u,C,W 2 ) 
Jo 

- [(pi+P2)-n]- N (7T(p 1 )7r(p 2 )\iljn- 1 (n-V) N ' 1 ^\0) 



A '- 1 " +J ,i/2 f n-(p 2 - Pl ) 

n n ; :(! K n-(p 2 + Pl ) 



X 

/o 



/ du6u(l-u)(2u-l) N - 1 Cl /2 (2u-l). (69) 

Jo 



The matrix elements of the local operators are related to each other by usual crossing 
symmetry, 

{p'\4)n--i(n- V)^" 1 ^) = (p,-p'\ipn--f(n-V) N - 1 i)\0). (70) 

[In this shorthand expression, both sides are regarded as functions of the pion four- 
momenta, defined in the respective physical regions, and analytic continuation is implied.] 
Using this relation with Eqs.(|68|) and (|69|) we obtain Eq.(|67|). Note that on the R.H.S. 
of Eq.(|S7|) the coefficients B n i(t) are taken at negative argument (t < 0), whereas in the 
expansion of the two-pion DA, Eg . (|B~T1) , they are defined for positive W 2 . The corre- 
sponding analytic continuation can be accomplished with help of dispersion relations (see 
Ref. p[ for details), 

g^ k! at fc -k J Ami s {s — t — ^0) 

where 5f(s), (J = 0, 1) are the im scattering phase shifts in the isospin and 1 channels. 
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Let us see the implications of the crossing relation, Eq. fl6"T|) , for the lowest moments 
of the skewed distribution. For the first moment of the isovector distribution we have 

J^dXH^iX^t) = B 01 (t) = F^it), (72) 

where F%' m '(t) is the pion electromagnetic form factor in the spacelike region, in agreement 
with Eq.fllB"!). For the second moment of the isoscalar distribution we obtain, substituting 
the explicit form of the Gegenbauer polynomials in £/2: 



f dXX H I=0 (X,^t) 



-1 Q 

5 



f 2 12 - £ 2 



(73) 



At t — 0, using the soft-pion theorem, Eq . (|64D , which implies B 1Q (0) + B 12 (0) = 0, we 
get 

£ i dXXH I=o (X 1 t;t = 0) = ^12(0) (l-j)- (74) 

If we substitute the value -812(0) = 5/9, which was computed in Ref. ||T9"f within the in- 
stanton vacuum model, we obtain precisely the generalized momentum sum rule, Eq . (|12"D . 

4.3 Application to quark/antiquark distributions in pion 

An particularly interesting application of the crossing relation, Eq.([H]), is the forward 
limit, t = and £ = 0, where the skewed distributions reduce to the usual quark/antiquark 
distributions, cf. Eq . fllPf) . The quark/antiquark distribution are known with fair accuracy 



from parametrizations of 7tN Drell-Yan and other data [[31]]. Eq.fl67|) relates the moments 
of the quark/antiquark distributions in the pion to the expansion coefficients B n \ at t = 0, 



which can in principle be measured in two-pion production at low invariant masses |T9 
The relation takes the form [cf. Eq. ([T0|) 1 

,1 V1 f <hW } f B&-W0) N odd, 

M? = / dXX N - 1 \ = A N { (75) 

1 i?r i,iv(0) A even, 




where the A^ are numerical coefficients which can be determined from Eq . (|67|) : A\ = 
1, A 2 = 9/5, A3 = 6/7, A 4 = 5/3, etc. For the lowest moments Eq. (ff5|) implies 
5q^ =1 (0) = M[ n = 1, which corresponds to the normalization condition Eq. (|59"D , and 
B{=° (0) = 5/9 Mf } , which corresponds to Eq.(|rj).g 

A non-trivial relation is obtained for N = 3. Using Eq.([75]) and the soft-pion theorem 
for the isovector two-pion DA, Eq. (|63|) , we can determine the coefficient -B2i(0) describing 
the deviation of the isovector two-pion DA from its asymptotic form, cf. Eq. (|62"D . One 
finds 

Bif(0) = a? - ? -Mt\ (76) 



3 To see this one needs to use the soft-pion theorem for the isoscalar two-pion DA, Eq. (^ 
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where a 2 is the expansion coefficient for the pion DA, Eq.(|65|). This relation is model 
independent an can be used as a consistency check for model calculations. 

Computation of the valence quark distribution in the pion in the effective theory based 
on the instanton vacuum (see Section |3|) gives a value of M3 = 0.25. [In this calculation 
the instanton-induced form factors F(d 2 ) have been taken into account.] The second 
moment of the pion DA has been computed in the same approach in Ref.|T6|, a 2 = 0.062; 



this small value is consistent with the CLEO measurements [pop . Substituting these results 
in Eq.(|76|) we obtain B^ 1 ^) = —0.22, which agrees with the result of a direct calculation 
in the instanton vacuum in Ref.|l9|.[] Thus, we see that the results obtained from the 
effective theory based on the instanton vacuum are consistent with the relation Eq. (|76]) . 
This is an extremely non-trivial check, since it shows that this approach preserves the 
soft-pion theorems (i.e., chiral invariance) as well as crossing symmetry. 

We note that the value of M3 = 0.25 obtained from the instanton vacuum is some- 



what larger than that of the GRV parametrization at the low normalization point |31|| , 
Mg = 0.16, and in good agreement with the value extracted from QCD sum rules with 



rW 



non-local condensates by Belitsky, M3 « 0.29 [ p2jj . 

When the mr phase shifts in Eq . ([7T|) in the isovector channel are approximated by 
exchange of "elementary" resonances (p, p', . . .), it becomes possible to express the expan- 
sion coefficients of the isovector two-pion DA in terms of the moments of the distribution 
amplitudes of the resonances (see Ref.pi for details). Keeping only the dominant con- 



tribution from p exchange in Ref. [19 was obtained the relation 



a 



(p) _ R /=i/nU-CM, 



2 



Bif(0)e- CM ^ (77) 



where the coefficient C was estimated in the instanton vacuum, C ~ 0.6 GeV 2 . In this 
approximation Eq. (|76|) becomes 

a fe CM l = a^ - -M 3 W . (78) 

Apart from the value of C, which does not influence the sign of the L.H.S., this relation 
is again model independent, and we can use it as a test for models of resonance DA's. We 
already noted that the instanton vacuum predicts -B^f^O) ~ —0.2, and thus a negative 
value for a 2 , which is consistent with Eq. ([ft|) because of the small value for a 2 obtained 
in this approach; see above. This result for a 2 is in contradiction to the results of QCD 



sum rule calculations, both in the standard approach |33| and with non-local condensates 
3"3J| , which obtained positive values. Ref.|33| reported a value of a 2 = 0.18 ± 0.1. How- 



ever, this calculation appears to be consistent with Eq.(^g), since a comparable sum rule 
calculation of a 2 |35| arrives at a relatively large value of a 2 ~ 0.44, so that Eq.([ft|) is 



satisfied if one substitutes, say, the value for M% from the GRV parametrization |3T|, or 
a slightly larger one. On the other hand, QCD sum calculations with non-local conden- 
sates give a significantly smaller value for the second moment of the pion DA; Ref. [J36f 



9 The slight numerical difference with the result quoted in Ref. pi, B^i 1 ^) ~ —0.20, is due to different 
approximations used in treating the form factors F(d 2 ). 
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.M _ n n i e. c,™ „1p„ cf Ho 7 ^ d/fW 



estimates aJj = 0...0.15; see also Ref.||37||. M3 was estimated in a QCD sum rule 



calculation with non-local condensates in Ref. 32 , Mi fa 0.29. This could indicate that 



» 



<?> 



in the calculation of a 2 in Ref.|34]], which quotes a value of 0.08 ± 0.02, the error margin 



could be somewhat larger than estimated. 

5 Conclusions 

In this paper we have investigated the structure of non-forward matrix elements of light- 
ray operators at a low normalization point, and their representations in terms of skewed 
and double distributions. Our principal conclusion is that the skewed distribution gener- 
ally has a "two-component" structure, i.e., that one is dealing with essentially different 
functions in the "quark/antiquark distribution" region, X > £/2,X < — £/2, and the 
"meson distribution amplitude" region, — £/2 < X < £/2. This follows naturally from 
the double distribution representation of the matrix element if one explicitly includes all 
twist -2 structures in the double distribution representation. We have shown that the con- 
tributions resulting from the "longitudinal" (oc r M ) twist-2 structure cannot be obtained 
from delta function terms in the conventional double distribution, which have previously 
been discussed by Radyushkin || |6 . 



This conclusion concerning the behavior of the skewed distribution at a low normaliza- 
tion point does not depend on any assumptions about the details of the non-perturbative 
at low scales. In fact, we have found qualitatively similar behavior in two different dynam- 
ical models: i) the low-energy effective theory based on the instanton vacuum, and ii) a 
generic meson exchange model. Moreover, the necessity to include terms ex r M in the de- 
composition of the non-forward matrix element Eq.(|T9"D is revealed by crossing invariance, 
which relates these terms to the matrix element for production of two pions. 

Our results concerning the general structure of non-forward matrix elements can 
serve as a basis for the construction of realistic parametrizations of skewed distribu- 
tions, satisfying all known requirements, and reproducing the phenomenologically known 
quark/antiquark and gluon distributions as well as the form factors of local operators in 
the appropriate limits. We plan to address this topic in a future publication. Whether or 
not the double distribution representation, in its complete form, Eg. (|24]) , will prove to be 
a useful tool for modeling skewed distributions remains to seen. 

The crossing relations between the skewed quark distribution and the two-pion dis- 
tribution amplitude derived here can be used to obtain additional information about the 
quark/antiquark distribution in the pion from measurements of electroproduction of two 
pions in 7*7 and 7*iV reactions. This fundamental characteristic of the pion is up to know 
only poorly known, since it can be measured directly only in hadronic ttN reactions such 
as Drell-Yan production, where it enters always together with the (anti— ) quark distri- 
butions in the nucleon. In particular, two-pion production in 7*7 reactions provides an 
opportunity to measure the distributions in the pion in a purely electromagnetic process. 
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A Generalized momentum sum rule for the pion 

In this appendix we give a derivation of the generalized momentum sum rule for the 
isoscalar skewed distribution in the pion, Eq.flT^D, using the universal chiral Lagrangian. 
The sum of the second moments of the skewed quark and gluon distributions is related 
to the form factors of the energy-momentum tensor ||. The symmetric version of the 
QCD EM-tensor is given by 

7^(0) = ^(0) 7{(1 V^(0) + F^(0)F«(0) - \g,„F 2 . (A.l) 

The general form of its matrix element between pion states is 

(7r+(p-r/2)|T^(0) \7T + {p + r/2)) = A(t)p^p u + B{t){r„r v - r 2 g^ u ) (A.2) 

(t = r 2 ); all other structures vanish due to energy-momentum conservation, d^T^ v {x) = 0. 
Expanding the non-local operators in Eqs.(^) and (0) in the distance z, 

$(-z/2)z[-z/2,z/2]il;(z/2) = ^(0)^(0) + Vz>(0) 7 „ V„ V(0) + ... (A.3) 
z»z»F m {-z/2)[-z/2,z/2]F?(z/2) = z»z v F IJa (p)F?(0) + ... (A.4) 

and keeping in mind that r-z = C,p-z, one can easily show that 

J 1 _ i dXX[H I=0 + H G \{X^t) = l -A{t) + ^B{t). (A.5) 

At t — > the matrix element of the energy-momentum tensor can be computed from 
first principles, since soft pion dynamics is described by the universal chiral Lagrangian. 
From 

^soft-pion = t^tt) 2 (A. 6) 

one obtains 

Tj^soft-pkm = dfiirdvit - -g^u{d^) 2 . (A. 7) 

Calculation of the matrix element between pion states, cf. Eq. (|A.2j ), gives 

^soft- P ion(i = 0) = 2, -B soft _ pion (t = 0) = — -. (A. 8) 



Inserting this in Eq.( |A.5|) for t — > one obtains Eq.fll2l). 
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B Resonance exchange contributions to non-forward 
pion matrix element 

An interesting property of non-forward matrix elements of QCD operators is the possibility 
of "meson exchange" contributions. One may view the exchange of a t-channel resonance 
as a purely phenomenological model for the matrix elements at a low scale, in the spirit 
of the vector dominance picture for the pion and nucleon e.m. form factors |J38| . Also, 



at large-iV c , where QCD is believed to become equivalent to a theory of resonances, one 
may hope to eventually construct a complete description of the pion matrix element of 
QCD operators in terms of resonance exchange. 

In this appendix we compute the contribution of the i-channel exchange of a spin- J 
two-pion resonance to the non-forward matrix element of the twist-2 operator in the pion, 
Eq.(|l]). Specifically, we want to show that exchange of isoscalar resonances (even J) gives 
rise to the behavior of the matrix element stated in Eq.(|22|). 



The contribution to the pion matrix element of the light-ray operator, Eq.(19), from 
an exchange of a resonance of spin J is given by the amplitude for the pion to emit the 
resonance, the resonance propagator, and the matrix element for the resonance to be 
"absorbed" by the light-ray operator (see Fig.|2|). The coupling of the resonance to the 
pion is of the form 

(7r(p-r/2)R(r,J,X)\7c(p + r/2)) = g Rnir p^ . . .p pn e^* Pri , (B.l) 

with isospin structure analogous to that of Eq. ([]]). Here e^ v denotes the polarization 
tensor of the spin-J resonance, and the coupling constant, gR nn , can be related to the titi 
width of the resonance. The upper part of the diagram in Fig.|2| is the matrix element 

(0${-z/2)z i-z/2, z/2] i>(z/2)\R(r, J, A)) 
AX) ~v\ z vj r i 

= f R (M R ) J *"»' J'; / dy e*"' 2 Mv)- (B.2) 

yr-zy' i J-\ 

Here, <fij(y) is a twist-2 distribution amplitude of the spin-J resonance. Conservation of 
angular momentum implies that 

/>"«v)-0 for {°^| J -j/^: (B.3) 

Computing the diagram Fig.|2] we obtain for the resonance contribution to Eq.(|l]): 
M(p-z,r-z;t) R _ cxch = 



B 
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where Pj is the Legendre polynomial of degree J, which results from the contraction of 
the vectors p and z with the transverse projector 

p^...p^ Zvi ...z VJ Y J ^ )vi - vn ^:. Pn - ( R5 ) 

A 

In particular, in the limit pz — > the argument of the Legendre polynomial in Eq. ( |B.4j ) 
becomes zero. Since -Pj(O) = for odd J we see that the exchange of isovector (odd- J) 
resonances does not contribute to the value of the amplitude at p-z — > 0, in accordance 
with the property Eq.(|5]) of Ai I=1 discussed in Section |[ In the isoscalar case, however, 
Eq.( |B.4| ) is non-zero. This is precisely the behavior described in Eq. fl22|) , which makes a 
double distribution representation in the form Eq.(p^) impossible. In fact, in the modified 
representation of A4 I=0 (p-z, r-z; t), Eq. (pl|) , isoscalar exchange leads to a contribution 
described by 

nf , sr fR9R,AM R ) J (-1) J / 2 [(J - l)!!] 2 mfi , 

i^)*-exch = E M 2_ t ^73T)n ^(V)- (B.6) 

J even 

From Eq. ( |B.4|) one sees that exchange of a J = resonance ( "sigma meson" ) contributes 
only to the r-z-term in the double distribution representation, Eq. fl2J|) , and thus only to 
D(y), while even-spin resonances with J > 2 contribute to both F I=0 (x,y) and D(y). 
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Figure 1: The range of the variables x and y in the double distribution, Eq.(14). The 
reduction to the skewed distribution, H(X,£), is achieved by integrating the double dis- 
tribution over the line x + yt^/2 = X, cf. Eq . (|T6l) , shown here for the case that X > £/2 
(thick line). 
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Figure 2: Schematic representation of resonance exchange contributions to the non- 
forward matrix element of the twist-2 operator in the pion, Eq. ([!]). The upper blob 
denotes the distribution amplitude of the spin- J resonance, Eq.( |B.2| ). 
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(a) 



(b) 



Figure 3: Diagrams in the effective low-energy theory contributing to the skewed quark 
distribution at a low normalization point. The dashed line denotes the pion field, the 
solid line the quark propagator with the dynamical quark mass, [i$ — MF 2 (d 2 )]~ 1 , and 
the filled circles the quark-pion vertices contained in the effective action, Eqs.(f4l|) and 
(fi^), which include a form factor F(d 2 ) for each quark line. Diagram (a) contributes only 
to the isoscalar distribution, and vanishes in the forward limit (r — > 0). 
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Figure 4: The contributions from diagrams (a) and (b) (cf. Fig.^J) to the isoscalar skewed 
quark distribution in the pion, H I=0 (X,C,), at a low normalization point, as functions of 
X, for a value of £ = 1. (Here t = 0). Dashed lines: Results obtained neglecting the 
momentum dependence of the dynamical quark mass, cf. Eqs.(|50D and fl51"|). Solid lines: 
The corresponding contributions obtained when including the form factors, F(d 2 ). Note 
that the contribution from diagram (a) is non-zero only in the region — £/2 < X < £/2. 
The momentum dependence of the dynamical quark mass forces this contribution to 
vanish at the end points, X = ±£/2. 
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Figure 5: The total isoscalar skewed quark distribution in the pion, H I=0 (X, £), for £ = 1, 
being the sum of the two contributions (a) and (b) shown in Fig.|j. Dashed line: Result 
obtained neglecting the momentum dependence of the dynamical quark mass, cf. Eq.(|52|). 
Solid lines: Distribution obtained when including the form factors, F(d 2 ). Due to the 
vanishing of contribution (a) at X = ±£/2 the total distribution is continuous at these 
points. 
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